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We discuss a general definition of directional derivative of any tensor flow field and its practical 
applications in physics. It is shown that both Lagrangian and Eulerian descriptions as comple- 
mentary types of fiow field specifications adopted in modern theoretical hydrodynamics, imply two 
complementary types of directional derivatives as corresponding mathematical constructions. One 
of them is the Euler substantive derivative useful only in the context of initial Cauchy problem and 
the other, called here as the local directional derivative, arises only in the context of so-called final 
Cauchy problem. The choice between Lagrangian and Eulerian specifications is demonstrated to be 
equivalent to the choice between space-time with Euclidean and Minkowski metric for any flow field 
domain, respectively. Mathematical consideration is developed within the framework of congruen- 
cies for general 4-dimensional differentiable manifold. The analytical expression for local directional 
derivative is formulated in form of a theorem. Although the consideration is developed for one- 
component (scalar) flow fleld, it can be easily generalized to any tensor field. Some implications of 
the local directional derivative concept for the classical theory of fields are also explored. 



I. INTRODUCTION 

Physical background and motivation for Euler's math- 
ematical construction known as substantive (or material) 
derivative becomes manifest in the way how the whole 
flow field may be specified in hydrodynamics. Usually, 
two complementary types of specifications (or representa- 
tions) are thought to suffice in order to provide a general 
description of the flow field kinematics. The first, La- 
grangian specification is based on identifying individual 
elements or bits of fluid domain. This idea associates a 
fluid motion with a geometrical transformation Ht on the 
closure such that the set Htflg represents the same in- 
dividual bit of fluid at time t. This representation is valid 
only if the identiflcation can be maintained by some kind 
of labelling usually denoting the initial position at instant 
tfj. The second, Eulerian specification was conceived as 
dissociated from identification of individual bits of fluid, 
only making use of the flow quantities as functions of lo- 
cal position in space at each instant of time during the 
motion. This specification results especially useful for hy- 
drodynamics of liquids and electromagnetic field descrip- 
tion in which any attempt of Lagrangian identification is 
impossible. Thus, in Euler's approach the velocity vector 
field is a primary notion. Assigning a velocity vector to 
each point of the fluid domain, one obtains the system 
of ordinary differential equations and their solutions as 
integral curves, intimately related to the given velocity 
vector field. The Eulerian representation provides a time 
parameterization of the curve in local coordinate system 
as a differentiable mapping from an open set of into 

Both complementary types of specifications (Eulerian 
and Lagrangian), generally speaking, different and cer- 
tainly useful in complementary contexts, can be made 
mathematically equivalent under special conditions. In 
terms of modern notation, they result equivalent within 
the formulation of the initial Cauchy problem for an or- 



dinary differential equation for velocity field. It gives a 
kind of dictionary for translating from one specification 
to the other. 

In this respect it is interesting to note that from the 
very birth of theoretical hydrodynamics as an indepen- 
dent body of mathematical knowledge, the conventional 
formulation of the directional derivative of flow fleld 
quantities tacitly implies the equivalence with the La- 
grangian speciflcation related to the initial Cauchy prob- 
lem. Thus, according to this type of speciflcation, Euler's 
substantial (or material) derivative describes the rate 
of time variation of /-property of fluid element on its path 
from one to the other point of space. As Euler himself 
coined it two and a half centuries ago 1] , his mathemat- 
ical construction described the rate of time variation of 
material properties following the motion of the fluid. 

On the other hand, the conventional definition of the 
Lie derivative on general differentiable manifolds, also ad- 
mits the same interpretation related to the initial Cauchy 
problem within the framework of congruencies for a given 
parameterization {A^}. When comparing scalars, vectors 
or tensors at different points {A^} and {Xi + AA^} on a 
certain integral curve, those entities under comparison 
are Lie (or invariably) dragged back from {Ai -I- AAi} to 
the point {Xi}. This gives a unique difference and hence 
a unique derivative as the limit of the difference between 
the values of scalars, vectors or tensors at different points 
on a manifold. A notion of a Lie invariably dragged func- 
tion along a congruence Xi is essential in the conventional 
definition of the directional (or Lie's) derivative ^ along 
a flow (or general vector) fleld. It defines the rule neces- 
sary to compare values of a given mathematical entity at 
two different points. 

The differential operator ^ is a tangent vector to the 
curve r(Aj) on a manifold. This association of the concept 
of a directional derivative from the classical analysis al- 
lows to maintain a visual picture on the Lie derivative as 
a tangent vector that generates a kind of a motion along 
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the curve r(Ai). In fact, when a differentiable manifold 
is an ordinary EucHdean domain both Eulcr's and Lie's 
derivatives coincide. 

As we shall discuss in this paper, both Eulerian and 
Lagrangian complementary flow field specifications which 
are traditionally related to the initial Cauchy problem do 
not cover all possible situations. It motivates a definition 
of complementary mathematical object denoted here as a 
local directional derivative. As a counter-part of the con- 
ventional approach, it appears in this paper formulated 
within the so-called final Cauchy problem and describes 
the time variation of /-property at some fixed point in 
a local coordinate system. It allows the Euler descrip- 
tion of the flow fleld to be dissociated from any need of 
identification of individual bits of fluid and, hence, to 
be considered properly as function of a position at every 
instant of time. 



II. EULER'S DIRECTIONAL DERIVATIVE 
AND INITIAL CAUCHY PROBLEM 

As the first task, let us explore the relationship be- 
tween the Eulerian and Lagrangian specifications within 
the conventional definition of the directional derivative. 
To simplify our analysis, the further discussion will be 
based on the consideration of one-component (scalar) 
fluid moving in a three-dimensional Euclidean domain. 
We denote by /(r,t) some regular function in an ar- 
bitrary space-time coordinate system. For instance, in 
fluid dynamics or elasticity theory it could be a density 
of some physical medium p in Cartesian coordinates. 

If in Lagrangian specification a geometrical transfor- 
mation Ht represents a mapping of the closure J7o onto 
HfClo for the same individual bit of fluid at time t, then 
Ht also represents the function 0: 



specifications become mathematically equivalent when 
an initial condition Tq — r(tQ) for Cauchy problem is 
added to the ordinary diff'erential equation H2.3|l . It gives 
a rule for translating from one specification to the other. 
If any quantity has the Eulerian representation /(r, i), 
its Lagrangian representation is : 



9(ro,t) = /(r(ro,t),i) 



(2.4) 



and, therefore, Euler's material or substantive derivative 
is conventionally defined as 0: 



Dt 



d 



d 



(l^ + (V'V))/ = -5(ro,i) 



df 



(2.5) 



where Tq is fixed. 

The differential operator has meaning only when 
applied to flow field variables as functions of (r(i), t) and 
gives the definition of directional derivative as a time 
derivative following the motion of the fluid in the di- 
rection of its velocity field v: 



D 
Dt 



/(r(t),t) 



1 



lhn-[/(r(to + i),to + t) 



/(r(io),to)] 
(2.6) 

Here, according to the standard definition H2.5|l related 
to the initial Cauchy problem r(<o) = i"o, both values 
/o = f{r{to),to) and / = f{r{to + t),to + t) represent 
the /-property at two different points of space Tq and 
ro -I- Ar, respectively. 

Thus, in this type of Eulerian of specification, one can 
reconstruct the property / of any identified bit of fluid 
at a new position r(io + dt) = v^ + dr and instant t^ -\- dt, 
based only on the knowledge of partial time derivative ^ 
and local distributions of the gradient V/ and velocity 
field V in local coordinate system: 



r = HtTa =^ r(rQ,t) 



(2.1) 



where points r and Tq denote the position- vector of the 
fluid identifiable point-particle at time t and initial time 
to, respectively. Thus, the velocity vector field on the 
domain is defined as: 



v=^r(r„,t) 



(2.2) 



where Tq is fixed. 

In the context of Eulerian specification, there is no 
explicit consideration of the function r = r{r^,t). The 
primary notion is the velocity field 



dr 

dt 



= v(r,i) 



(2.3) 



as function of position in space (r) and time (t) on a fiuid 
domain. Based on this convention, it is clear that both 



/ = /o + (|j- + (v,V/)dt 



(2.7) 



The function / has, generally speaking, explicit as well 
as implicit (through r(t)) time dependencies and, there- 
fore, may be defined on a 4-dimensional space-time man- 
ifold with no metric known a priori. In this respect, 
the Lie derivative, as a particularly interesting general- 
ization of H2.5|l (or (|2.6() '). will be convenient for further 
considerations, since it provides a necessary framework 
on manifolds without metric. The differential equation 
(|2.3|) will define a congruence or i-parameterized set of 
integral world-lines filling a 4-dimensional manifold: 



dx' 
'dt 



to+t 

V'; x\t) =4+ j V'dt 

to 



(2., 



where x = (a;°, a;^, a;^, a;"^) = (t,r); V = (l,v); Xq = x{to) 
and, for our convenience, we leave for the time variable x^ 



3 



its original denomination t. Upper indices are used for 
the coordinate functions a;*(t) so that the 1-forms will 
satisfy the index conventions of modern differential ge- 
ometry. 

If the velocity vector field V is C°°, the coordinate 
transformation H2.8|l is a diffeomorphism, forming part 
of a one-parameter Lie group. Let us denote this trans- 
formation as Ft: (xq —>■ x{xo-,t)), which defines the 
mapping of f{xo) = /(io,r(^n]) along the congruence 
(called also as Lie dragging Q) into a new function 
f{xo + x) = f{to + t,r{to + t)). A Lie dragging of 
scalar field has a simple geometrical interpretation in La- 
grangian specification of the fluid field: Ft transforms the 
/-property of the identified fluid element at Xq according 
the rule j^]: 



iFtf)ixo)^ fiFtixo))^ fixo + x) 



(2.9) 



into the /-property of the same fluid element at xq + x. 
This interpretation also concerns the analytic expression 
of the Lie derivative Ly along the vector field V ~ (1, v): 



Lvf = [^Ft .f]t„ = lim i[/(.To + x) - f{x^)] (2.10) 
at t^o t 

where x{tQ) = 0. The concept of a Lie invariably dragged 
function along a congruence is used in this conventional 
definition (and implicitly in (|2.6|l ): in fact, the quantity 
f{xQ+x) is invariably dragged along the congruence from 
xq-\-x back to xq, since in general tensor calculus it has no 
clear meaning to compare both values f{xo) and f{xo+x) 
at different points in a manifold without metric. Thus, 
(|2.10() gives a unique difference and therefore a unique 
derivative. In fact, when t is too small, the mapping Ff. 
{xq — > x{x(),t)) has an explicit form: 



Ftxo = x'{xo,t) = xl + tV^xo) + o{t) (2.11) 
which gives analytic expression for the Lie derivative: 



76 (2^12) 



Important that both traditional definitions for Euler's 
and Lie's directional derivatives, respectively, turn out 
to be defined entirely in the spirit of original Lagrangian 
specification, i.e. when a fluid element or a point on a 
congruence are constantly identified in a local coordinate 
system. It explains why in an ordinary Euclidean domain 
the Lie mathematical construction takes a familiar form 
of Euler's directional derivative 0: 



V)/ = ^ 
Dt 



(2.13) 



It is known as a full derivative along the vector filed 
V = (l,v) and in hydrodynamics it also has numerous 



applications for the description of the motion of macro- 
scopic individual bodies of fluid. 

In this respect, let us now consider a time variation of 
the fluid /-content in a macroscopic 3-dimensional space 
domain V moving with the fluid: 

iff^y-ii: I f^y-EiuSV) (2.14, 



V{t) 



SV{t) 



where, for our convenience, the macroscopic volume V{t) 
is represented as the sum of individual microscopic vol- 
umes SV{t), i.e. V = J2SV. 

A geometrical transformation Ht describes the evolu- 
tion of the whole volume V{t) = HtVo as well as its 
individual bits SV{t) = Ht6Vo. The assumption that the 
domain V and all 6V move with the fluid, means that 
there is no flux of / through the common bounding sur- 
face of the fluid domain dV and bounding surfaces of all 
individual elements dSV. It imposes the condition that 
normal components of the relative fluid velocity field are 
zero at every point of all such surfaces moving with the 
fluid. By the chain rule we obtain from H2.14|) : 



V(t) 



V[t) 



The time variation of a microscopic volume 5V is a 
result of movement of each point of the bounding surface 
dSV and it is described by the divergence theorem Q: 



j^5V{t)= J {n,u)dS = J{V,u)dV 



(2.16) 



dSV 



SV 



where n is a unit vector normal to the surface dSV and 
u is the local velocity field of surface points in the co- 
ordinate system attached to SV. Since the divergence 
does not depend on the choice of coordinate basis, we 
can write (V,u) = (V,v) for the flow velocity field v 
defined in the main coordinate system at rest. Hence 
the local rate of expansion or dilation of a microscopic 
volume element 6V is: 



-SV{t)^iV,^r)^V 



(2.17) 



When this is substituted into (|2.15|) . one gets the result 
of the Convection Theorem: 



d_ 
'dt 



fdV 



Dt 



/(V,v))dl/ 



(2.18) 



V(t) 



v(t) 



More detailed and rigorous demonstrations can be found, 
for instance, in 0, We only need the way it was rea- 
soned as well as its interpretations for further discussions. 
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To finish this Section, we conclude that the formula 
(2.13) for Euler's derivative (as well as its integral 
counter-part H2.18|l l provide the rule for reconstruction 
of the /-property (or the /-content) at a new position 
r = ro -f dr (or in a new domain V{t) = HtVo), gener- 
ally speaking, different from the initial one. This hydro- 
dynamics interpretation will become useful in the next 
Section in order to contrast a complementary specifica- 
tion of a flow field which provides the reconstruction of 
the /-property at a fixed point of space Vq but at instant 
to + dt different from tg, entirely in the spirit of Euler's 
original idea of flow field specification. 

III. FINAL CAUCHY PROBLEM AND LOCAL 
DIRECTIONAL DERIVATIVE 

Let us discuss in this Section some limitations of con- 
ventional definitions (|2.6|) and (|2.1UI) for a flow field in 
situations when Lagrangian identification is in principle 
impossible. Before proceed, it will be convenient to ap- 
ply the above-used terminology to clarify the definition 
of the partial time derivative of a fiow field /-quantity in 
the framework of the classical analysis: 

§f = liml[/(ro,to + t)-/(ro,io)] (3.1) 
dt t-»o t 

where note that ^ means not only the rate of change of 
/ at a fixed position ro in space but (|3.1|) is considered 
under the condition that the space variable ro is not a 
subject of time parameter t at all (i.e. ro(t) = const). 
In other words, this situation corresponds to the time 
derivative of flow quantities with a frozen velocity field 
V (i.e. V = 0). 

Let us now consider a time variation of the fluid /- 
content in a fixed 3-dimensional space domain Vq. If we 
follow the way of reasoning used earlier for the formula- 
tion of the Convection Theorem, then: 




where again, for our convenience, the macroscopic vol- 
ume Vq is represented as the sum of individual micro- 
scopic volumes SVq, i.e. Vq = J2 ^^o- 

Since all volume elements SYq are now fixed, the differ- 
ential operator ^ acts only on /-property and is used to 
be associated with the partial time derivative ^, giving 
place to the well-known relationship for a fixed domain 
Vq (see any text on the classical theory of fields): 




Vo Vo 

A note of caution is appropriate here. One might sus- 
pect that some type of specification of the flow fleld for 



/-property inside SVq is necessary as it was in (I2.15|) for 
the reliable use of j^. The integrand in the left-hand 
side of H3.2() is not a simple multivariable function /(r, t) 
but a rather different mathematical entity f{r{t),t), i.e. 
flow field quantity. This feature makes it difficult the 
straightforward application of the partial time derivative 
^ meaningful for a fixed space variable r as function of 
a frozen velocity field (v — 0). Therefore, the use of 
the partial time derivative does not seem to be fully jus- 
tified in the conventional approach. This circumstance 
was also critically pointed out in 0-0 ■ 

To clarify the situation, let us go through the same 
mathematical construction as was used before for Euler's 
and Lie's derivatives, i.e we shall define an appropriate 
specification of the flow field in this case. The terminol- 
ogy of ordinary differential equations theory, combined 
with notions from classical analysis, enables us to give 
a useful and compact definition of the left-hand side of 

Let ro be a fixed point of the closure Vq and let the 
path of some identified elementary bit of fluid (that at 
some earlier instant to passed through a certain point of 
space r*{to)) lie at present instant <o + ^ on the position 
of the fixed point of space ro . Then the full time deriva- 
tive is understood as the limit of the difference between 
the values of the volume /-content at different instants 
to and to + t. This requirement provides a framework 
necessary to derive analytic expression for the left-hand 
side of l|X^ : 

^ J f{v{t),t)dV^ 

Vo 

lim i j [f[v{to + t), to + t)- /(r(to), to)] dV (3.4) 

Vo 

where r(to) = ro. Both values /(r(to),to) and /(r(to + 
t) , to + t) represent the /-property at the same point of 
space ro. 

Based on this convention we note that r(to + t) = ro 
should lie at the end of the integral curve: 

to+t 

ro = r*(to)+ j ^^{v,t)dt (3.5) 

to 

Extrapolation on a set of integral curves that fill our 
domain is straightforward. The path of every such curve 
has to end at some fixed point of the closure Vq. All 
of them are solutions of initial Cauchy problems for the 
first-order differential equation: 

dr 

- = v; r(to) = r* e (3.6) 

In Lagrangian specification it represents a geometri- 
cal transformation (or mapping) Ht of the original 3- 
dimensional domain V* at instant to onto HfV* = Vq at 
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instant to+t. If the inverse mapping HjT^ is single- valued, 
then V* = H^^Vq can be regarded as a reconstruction of 
initial conditions from the knowledge of the final domain 
Vo, i.e. V* becomes dependent on t parameter. Simi- 
larly to (|2.1|l . for an individual bit of fluid, Ht defines 
the function: 



The inconvenience of this description is that we are now 
at the local coordinate system as a function of r*(t). 
Since the integration in 1)3. 4|l is effected over the fixed 
domain Vb, we choose a local coordinate system attached 
to Vb by means of coordinate transformation: 



to+t 



(3.7) 



where rg is fixed and belongs to Vq. 

By analogy with the initial Cauchy problem for H2.3|l . 
it may be called as a final Cauchy problem for H3.6|l . If 
any flow field quantity has the Eulerian representation 
/(r, t), its Lagrangian representation g* in the context of 
the final Cauchy problem will be: 



g*(ro,t) = /(r(ro,t),i) 



(3.8) 



and, therefore, the partial time derivative (with the 
fixed To) will define some new mathematical construction: 



9 ^ D*f 



(3.9) 



Note that the differential operator makes sense 
only when applied to flow field variables as functions of 
the final Cauchy problem and will be called here as the 
local directional derivative by analogy with the definition 



of -Sj within the framework of classical analysis: 



^=limi[/(r(to+t),to 



<)-/(r(to),to)] (3.10) 



where r(io) = Tq and we denote by r(to -|- t) the final 
point ro of the integral curve H3.5II . 



To find the analytic expression for , let us use a con- 
sideration similar to the applied in the previous Section 
for the definition of The differential equation (|3.6|) 
will define a congruence or ^-parameterized set of integral 
world-lines filling a 4-dimensional space-time manifold: 



dt 



to+t 



V'dt 



(3.11) 



r*(0=ro- J v{r,t)dt 



(3.13) 



or in 4-dimensional notations: 



dx* 



to+t 



dt 



= V* 



V"dt 



(3.14) 



to 



where V* = (1,— v). It defines an equivalent mapping 
which we denote as G^: (xq ^ x* = x{xQ,t)). 

Note that Gj is the same mapping Gt but defined in 
the local coordinate system attached to Vq. More pre- 
cisely, it means that the final Cauchy problem (|3.6II ad- 
mits an equivalent formulation as an initial Cauchy prob- 
lem (|3.13|) (i.e. Gj is not the inverse transformation 
(Gt)^'^, since the course of time is not changed on the 
opposite one). 

Thus, when t is too small, the transformation Gj has 
an explicit form: 



G;xo = x*\xa,t) = 4 + tV*\xo) + o{t) (3.15) 

which gives analytic expression for the Lie derivative 
along the congruence V* at a local coordinate system: 



/ 



dt 



f{G; xo) = V* 



5/ 

dx^ 



(3.16) 



In an ordinary Euclidean domain this mathematical con- 
struction takes the following expression: 



V^/ = (|-(v,V))/ = f;{ (3.17) 

Now, noting that the integrand in the right-hand side 
of the equation (|3.4() is jpj according to the definition 
(|3.1U|) . we can proceed to the formulation of our result 
(|3.17|) as a theorem proven in an ordinary Euclidean do- 



where X* = {to, r*) and xq = {to + t,ro). 

This transformation, which we denote as Gt'. {x* —> 
Xq = x{x*,t)), defines the mapping of f{x*) along the 
congruence into a new function f{xo). In Lagrangian 
specification of the fiow field it has an obvious geometri- 
cal interpretation: Gt transforms the / property of the 
identified fluid element at x* according the rule: 



{Gtf){x*) = f{Gt{x*))=f{xo) 



(3.12) 



Theorem 1 (Local Convection Theorem): Let w he a 
vector field generating a fluid flow through a fixed 3- 
dimensional domain Vq and if f{r,t) G C^{Vo), then 



d 
dt 



J ^"^^^ y'(^-(v,v))/dF 



(3.18) 



Vo Vo 
where dV denotes the fixed volume element. 
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This result formulated within the framework of final 
Cauchy problems could be regarded as a complementary 
counter-part of the Convection Theorem H2.18|l consid- 
ered within the framework of initial Cauchy problems. 

In Eulerian specification, the formula (|3.17|) admits a 
clear hydrodynamic interpretation: provides the rule 
for reconstruction of the /-property at a fixed point of 
space rg at instant to + dt, based only on the knowledge 
of partial time derivative ^ and local distributions of 
gradient V/ and velocity field v in the vicinity of Tq: 



Dt 



D*f 
D*t 



(3.23) 



where gij = Sij is the Euclidean metric tensor. A 
Minkowski metric is also consistently singled out for local 
directional derivative 



D*f 

D*t 



(3.24) 



/ = /o + ( 



a/ 

dt 



Vf)dt 



(3.19) 



Generally speaking, this type of Eulerian specification 
of the flow field does not imply any sort of identification 
of fluid elements and hence ought to be complementary 
to the original Lagrangian approach. In fact, it considers 
the rate of time variation of /-property locally, at fixed 
position of space. Whereas the Euler derivative comple- 
mentarily describes the rate of time variation following 
the motion of the fluid. 

Both types of directional derivatives and -^^j can 
be analyzed in terms of 1-forms or real-valued functions 
of vectors in 4-dimcnsional manifolds: 



(3.20) 



where V* = (1, — v); g*^ = diag{l, —1, —1, —1) is indefi- 
nite or Minkowski metric tensor. 

Another consequence of this form is that it gives or- 
thonormal bases for space-time manifolds (previously 
introduced with no metric known a priori). For La- 
grangian flow field specification, a basis is Cartesian and 
a transformation matrix A^ from one such basis to an- 
other is orthogonal matrix: 



A = A: 



9i3 



= A: 



9ij 



A, 



(3.25) 



These matrices Ac form the symmetry group 0(4). 

Likewise, for Eulerian specification a Minkowski met- 
ric picks out a preferred set of bases known as pseudo- 
Euclidean or Lorentz bases. A transformation matrix Al 
from one Lorentz basis to another satisfies: 



where « = 0, 1, 2, 3 and (^r) = (f^, V/) in an ordinary 
Euclidean domain. 

Now we point out that in tensor algebra the set {ujiV^} 

are components of a linear operator or ( | ) tensor. The 

formation of a scalar oj{V) is called the contraction of 
the 1-form uj with the vector V and it is an alternative 
representation of directional derivatives: 



£1 

Dt 



D*t 



(3.21) 



The contraction of diagonal components of the tensor 
ujiV^ is independent of the basis. Importantly, this law 
shows that both types of directional derivatives and 

are invariant and do not depend on the choice of a 
local coordinate system. On the other hand, it is also the 
property of scalar product in manifolds with metric. The 
metric tensor maps 1-forms into vectors in a 1-1 manner. 
This pairing is usually written as: 



gij^ 



V' = g'Wj 



(3.22) 



Therefore, from the point of view of tensor algebra, (|3.21|) 
can be considered as a scalar product in a 4-dimensional 
manifold with metric: 



A^ 



5r 



A; 



(3.26) 



Ai is called a Lorentz transformation and belongs to the 
Lorentz group L(4) or 0(3, 1). 

The point that needs to be emphasized here is the re- 
markable circumstance of Euler's specification in evoking 
of the Minkowski metric without any previous postula- 
tion. In other words, consistent mathematical descrip- 
tion of fluids is perfectly compatible with the Lorentz 
symmetry group. This fact was not seriously considered 
in theoretical hydrodynamics until now. 

The Gahlean group as one of subgroups of 0(4), is 
commonly used in modern classical mechanics in flat 
space-time manifolds. This is not surprising in view that 
all classical mechanics laws are written in Lagrangian 
specification by constant identification of mechanical ob- 
jects and within the formulation of initial Cauchy prob- 
lem for equations of motion. It was therefore natural 
to admit that space-time in classical mechanics has a 
Galilean group symmetry. The Special Relativity postu- 
lation of Lorentz group symmetry on mechanics is not 
trivial, having in mind the complementary character of 
Lagrangian and Euler's descriptions. Perhaps it can ex- 
plain a paradoxical nature of some conclusions in rela- 
tivistic mechanics but it overcomes the scope of this work 
and will be considered elsewhere. 

Thus, from the point of view of flow field specifications, 
both kinds of directional derivatives are complementary 
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and equally valid but should be used in different contexts. 
Euler's derivative has therefore a more narrow framework 
of applicability in the classical field theory than it was 
supposed. In what follows we will confine our attention 
on some example from the classical field theory. 



IV. LOCAL DIRECTIONAL DERIVATIVE IN 
CLASSICAL FIELD THEORY 

Let us now consider the description of the conserva- 
tion of the fluid /-content in an arbitrary 3-dimensional 
space domain. If the volume V moves with the fluid, the 
Convection Theorem (2.18) written in differential form 
states that the /-content is conserved when the total time 
derivative is zero: 



Dt 



/(V,v) 



9/ 
dt 



V(/v) - 



(4.1) 



In particular, when the velocity field v is locally zero, it 
represents the continuity equation of any elastic medium 
locally at rest: 



dl 

dt 



-f /(V,v) = 



(4.2) 



where the extra term (v, V/) due to the fluid movement 
has disappeared. 

In the case when the fluid moves through a volume Vq 
fixed in local coordinate system at rest, a mathematical 
restriction on conservation immediately leads to the well- 
established integro-differential form of continuity equa- 
tion: 



)dV (4.3) 



^ J fdV j f{wAS) = - j V(/v 

Vo So Vo 



Note that here both sides of the equation are obviously 
independent on the choice of a particular coordinate ba- 
sis. 

As we already mentioned earlier in the previous Sec- 
tion, it is commonly thought that, in this case, the total 
time derivative can be substituted in the integrand by the 
partial derivative, giving place to the conventional form 
of continuity equation in the reference system at rest: 



(4.4) 



The circumstance that it coincides with the expression 
(|4.1|) derived for the volume in motion, is mainly at- 
tributed to the cross-verification of the standard dif- 
ferential form of continuity equation. Nevertheless, it 
is strange to contemplate that the differential equation 
(|4.4|) does not possess the symmetry properties of its 



original integral counter-part (|4.3(l . The left-hand side 
of (|4.4|) becomes manifestly dependent on the choice of 
a coordinate basis which, generally speaking, leads to a 
more narrow group of symmetries. Let us see whether 
results of the previous Section may help to clarify the 
situation. 

In fact, implementation of Euler's type of flow field 
specification for the left-hand side of (|4.3|l in the frame- 
work of the final Cauchy problem changes the character 
of the integrand expression. If it is considered as the 
local directional derivative of /-property, the continuity 
equation H4.3|l takes the following differential form: 



D*f df 



dt 



D*t 



(v,V)/ = -V(/v) 



(4.5) 



that coincides with (|4.2|) . The right-hand side of 1)4. 5|) 
as a divergence and the left-hand side as the local direc- 
tional derivative do not dependent on the choice of 
a coordinate basis. It means that this differential form 
of the continuity equation has the symmetry properties 
of its original integral counter-part (|4.3|l . 

A brief comment is worthy in this respect. Why the 
traditional approach based on the Convection Theorem 
gives a different result 1)4. l|l ? Certainly, it is correct 
but it has a non-invariant extra term (v, V/) due to 
the fact that the description is effected in the reference 
system at rest for the domain following the motion of 
the fluid. These shortcomings of the direct application 
of the Convection Theorem was not appreciated until 
now. If an observer moves with the fluid, this Theorem 
gives the equation H4.2|l without an extra term and with 
the symmetry properties of the original integral equation 
(14. 3|) . On the other hand, one could logically ask why 
all numerical simulations based on the standard differen- 
tial form of the continuity equation (|4.4|) do not lead to 
incorrect predictions? The answer is the following: tra- 
ditional time discretization schemes for the partial time 
derivative of flow field quantities (see, for instance, 0) 
treat it as if it were the total time derivative. 

Another interesting task would be an application of the 
local derivative concept to the integral form of Maxwell's 
equations. Two of them contain the full time derivative 
over volume integrals and are known as induction laws 
for electric E and magnetic H vector fields, respectively, 
in the local frame of reference at rest: 



(H,dl) 



An 
c 



(j,dS) 



1 d 
c dt 



{V,E)dV (4.6) 



(E,dl) 



1 d 

c dt 



{V,Il)dV 



(4.7) 



Straightforward application of H3.18|l in this case is hin- 
dered by a priori unknown nature of the velocity vector 
field for electric and magnetic field components. At this 
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stage, only quasistatic approximation can admit a reli- 
able application of the local directional derivative con- 
cept. In fact, the Special Relativity firmly established 
that electromagnetic field components of uniformly mov- 
ing single charge do not depend explicitly on time pa- 
rameter t. In other words, E and H are thought to be 
rigidly attached to the charged particle and uniformly 
move with it. This is one of the consequences of the 
Relativity Principle. Thus, if the charge velocity Vg is 
known, the velocity vector field v for quasistatic compo- 
nents of electric and magnetic field is also defined in the 
closure V. Applying the result of the Theorem H3.18|l . we 
can rewrite (|4.6I) - (I4.7II in a more convenient form: 

J (H, dl)^^J (j, dS) i y'(^ - (v, V))VEd^ 

C S V 

(4.8) 

I (E, dl) ^ -1 y (^ - (v, \/))VIldV (4.9) 

c V 

where v = Vq is the instantaneous velocity field in the 
closure V. 

Since the motion is uniform, all partial derivatives van- 
ish from (|4.8|) - (|4.9|) . Applying a well-known expression 
for a general vector field A: 

(v,V)A = v(VA)-[V,[v,A]] (4.10) 

and reducing the volume V to zero, we arrive to the 
well-established relationship between quasistatic mag- 
netic and electric field strength of an uniformly moving 
charge from the point of view of a local inertial reference 
system 

H = i[v,E]; E=-i[v,H] (4.11) 

c c 

It is worth stressing that a priori no relativity prin- 
ciple was needed in deriving these transformation rules 
for electric and magnetic field components. The term 
proportional to — (v, V) can be considered as convective 
displacement current [T^ by analogy with Maxwell's dis- 
placement current proportional to N ote that the in- 
tegral form of Maxwell's equations l|4.8|l - 14.9|) written in 
Euler's specification is now compatible with the charge 
conservation law (|4.5|) also represented in Euler's speci- 
fication. 

On the other hand, the Lorentz and Ampere force con- 
ceptions are manifestly valid quasistatic approximations 
(|4.11l) and therefore are inclosed into integral form of 
Maxwell's field equations on a basic level. It means that 
there may be no need to postulate them separately as it 
was done in Maxwell-Lorentz microscopic electron the- 
ory and remains accepted at present. Nevertheless, any 
full analysis of these issues comes out of the scope of the 
present consideration and will be given elsewhere. 



V. CONCLUSIONS 

We attempted to consider a logical background and 
structure of useful mathematical constructions which are 
traditionally based on both Eulerian and Lagrangian flow 
field representations, complementary to each other. This 
account provides a mathematical method that justifies 
the definition of a complementary counter-part for Eu- 
ler's directional derivative which is called here as the local 
directional derivative. 

The point that needs to be emphasized is the comple- 
mentary character of the above introduced concept. By 
no means it substitutes the Euler mathematical construc- 
tion. By contrary, it is shown that both types of direc- 
tional derivatives are equally valid but should be used 
in different contexts. In fact, Euler's substantive deriva- 
tive arises in the context of initial Cauchy problems and 
therefore becomes useful within the framework of the La- 
grangian type of description of flow field quantities. Like- 
wise, it is possible to define a complementary framework 
of so-called final Cauchy problems appropriate for the 
Euler flow field specification as a function of position in 
space and in time for fiuid domain. From the point of 
view of the classical theory of fields it means a more nar- 
row framework of applicability for Euler's derivative than 
it was thought. 

The analytic expression for the local directional deriva- 
tive is formulated in form of a theorem analogous 
and complementary to the Convection Theorem well- 
established in theoretical hydrodynamics. One of its 
interesting conclusions is that the choice between La- 
grangian and Eulerian types of fiow field representation 
is equivalent to the choice between space-time mani- 
folds with Euclidean and Minkowski metric, respectively. 
Therefore, the consistent mathematical description of 
kinematics of fluids in Eulerian representation results 
compatible with the Lorentz group symmetry L{4). In 
fact, it could be understood as complementary to the 
0(4) group symmetry compatible with the Lagrangian 
representation. 

On the other hand, the definition of the mathematical 
construction complementary to the traditional one, helps 
to get a deeper insight on the cross- verification of several 
partial differential equations obtained from their well- 
established integral counter-parts in classical theory of 
fields. 

Although the consideration in this work was developed 
for one-component (scalar) fiow field, the notion of the 
local directional derivative can be easily generalized on 
Lie's derivatives for any general tensor field on differen- 
tiable manifolds. Both types of Lie's derivative will cor- 
respond to both complementary types of specifications. 

In place of concluding remark let us give asserting and 
encouraging words of a great mathematician. Gauss once 
wrote in his letter to Bessel (quoted form ---One 
should never forget that the function [of complex 
variable], like all mathematical constructions, are only 
our own creations, and that when the definition with 
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which one begins ceases to make sense, one should not ticular, to Prof. A.E. Chubykalo, Prof. V. Kassandrov 
ask, what is, but what is convenient to assume in order and Dr. V. Onoochin. 
that it remain significant..." . 
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